arXiv:1501.04356v2 [cond-mat.str-el] 10 Feb 2016 


S.-Z. Lin and A. Saxena , Phys. Rev. B 93, 060401(R) (2016). 


Dynamics of Dirac strings and monopole-like excitations in chiral magnets under a current drive 

Shi-Zeng Lin* and Avadh Saxena 

Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 

(Dated: February 11, 2016) 

Skyrmion lines in metallic chiral magnets carry an emergent magnetic field experienced by the conduction 
electrons. The infiow and outfiow of this field across a closed surface is not necessarily equal, thus it allows 
for the existence of emergent monopoles. One example is a segment of skyrmion line inside a crystal, where 
a monopole and antimonopole pair is connected by the emergent magnetic flux line. This is a realization of 
Dirac string-like excitations. Here we study the dynamics of monopoles in chiral magnets under an electric 
current. We show that in the process of creation of skyrmion lines, skyrmion line segments are first created via 
the proliferation of monopoles and antimonopoles. Then these line segments join and span the whole system 
through the annihilation of monopoles. The skyrmion lines are destroyed via the proliferation of monopoles and 
antimonopoles at high currents, resulting in a chiral liquid phase. We also propose to create the monopoles in 
a controlled way by applying an inhomogenous current to a crystal. Remarkably, an electric field component 
in the magnetic field direction proportional to the current squared in the low current region is induced by the 
motion of distorted skyrmion lines, in addition to the Hall and longitudinal voltage. The existence of monopoles 
can be inferred from transport or imaging measurements. 

PACS numbers: 75.10.Hk, 75.25.-j, 72.25.-b, 75.78.-n 


Skyrmion in magnets is a stable mesoscopic topological ex¬ 
citation where spins wrap the surface of a sphere once. [1,2] 
Skyrmions have been observed in magnetic materials, such as 
B20 chiral magnets without inversion symmetry [3-5], mul- 
tiferroics [6, 7] and even centrosymmetric compounds [8, 9]. 
These findings suggest that skyrmions are ubiquitous in mag¬ 
netic systems. The typical size of skyrmions is about 10 nm 
and skyrmions form a triangular lattice. Skyrmions can be 
driven by electric current [10-12], temperature gradient [13- 
15] and electric field gradient [16, 17]. Moreover the thresh¬ 
old current that makes skyrmions mobile is 5 to 6 orders of 
magnitude smaller than that for a magnetic domain wall. [ 10- 
12] Thus skyrmions are promising candidates for spintronics 
applications and recently they have attracted significant atten¬ 
tion. [18, 19] 

In metals, the spin of the conduction electron interacts with 
the localized magnetic moment associated with skyrmions 
through the Hund’s coupling. In the strong coupling limit, 
the spin of a conduction electron changes adiabatically when 
it moves around a skyrmion and it gains a Berry phase which 
can be described by an effective magnetic field. The emergent 
magnetic field is related to the skyrmion topological charge, 
= hcEijkn • (djU X dkn)l(2e), and is extremely strong due to 
the small size of the skyrmion. Here n is a unit vector repre¬ 
senting the direction of localized spin and eijk is the totally an¬ 
tisymmetric matrix with /, j, k = x, y, z. A skyrmion line in 
a 3D metal thus can be regarded as an emergent magnetic flux 
line with quantum flux Oq = hc/e. It was demonstrated ex¬ 
perimentally that two skyrmion lines can merge and become a 
single line. This implies an emergent magnetic monopole with 
flux Oo at the merging point, which can be verified by integrat¬ 
ing over a closed surface around the merging point. [20] 
A monopole and antimonopole pair can also be created by 
thermal fluctuations. [21] The emergent magnetic flux of the 
skyrmion and hence the monopole is an effective description 
for the conduction electrons [22], thus these concepts can¬ 


not be applied to insulators. While magnetic monopoles as 
elementary particles have never been found experimentally, 
quasiparticles resembling monopoles have been identified in 
several condensed-matter systems, such as spin ice [23] and 
chiral magnets discussed here as well as Bose-Einstein con¬ 
densates [24]. A unique feature of monopoles in skyrmion 
systems is that they can be driven by electric currents. The 
monopoles in chiral magnets have attracted considerable in¬ 
terest recently. [25-27] 

In a thin film, skyrmion becomes a pancake-like object and 
no monopole is allowed to appear because the system is uni¬ 
form along the direction perpendicular to the film. It was 
demonstrated theoretically that skyrmions can be created dy¬ 
namically by applying a strong current to a fully spin po¬ 
larized state in chiral magnetic thin Aims. [28] This is be¬ 
cause the current induces a Doppler shift in the magnon spec¬ 
trum [29, 30] and renders the spin polarized state unstable 
when the spectrum becomes gapless, after which skyrmions 
are nucleated dynamically. When the current increases fur¬ 
ther, skyrmions are destroyed and the system evolves into a 
chiral liquid with strong spatial and temporal fluctuations in 
spin chirality. 

The creation and destruction of skyrmions by current can 
also occur in 3D, where the emergent monopoles play an 
important role. A skyrmion line percolating the whole sys¬ 
tem (connecting surfaces of the system) cannot be created in¬ 
stantly. Instead segments of skyrmion lines are created. At the 
ends of these lines, V • B^ = +p^, i.e. there is a monopole and 
antimonopole pair connected by a skyrmion line, as schemat¬ 
ically shown in Fig. 1 (a). This is a realization of Dirac 
string in condensed matter systems. These segments then join 
together to form skyrmion lines spanning the whole system 
through the annihilation of monopoles and antimonopoles. In 
the course of skyrmion destruction, skyrmion lines break into 
segments and monopole and antimonopole pairs are created 
dynamically by current. These physical pictures are borne out 
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FIG. 1. (color online) (a) Schematic view of a skyrmion segment as a Dirac string-like excitation, where there is a monopole and antimonopole 
pair at the two ends, (b) Straight skyrmion line lattice, (c) Snapshot of spin configuration after the magnonic instability of spin polarized 
state. Skyrmion segments are created dynamically associated with monopole and antimonopole excitations. The monopole and antimonopole 
annihilate or disappear at the surfaces, through the merging of the skyrmion segments. Finally the system reaches a state with straight skyrmion 
line lattice, and they move as a whole, (d) Chiral liquid at a high current when skyrmion lines are destroyed. Color at three surfaces represents 
Uz (red for n^- I and blue for = -1) and the vector field at the top surface denotes the and Uy components. 


by our numerical simulations below. 

We consider a 3D chiral magnet described by the following 
Hamiltonian in the continuum limit because the skyrmion size 
is much bigger than the underlying lattice parameter [ 1 , 2 , 31- 
33] 


7Y 


-M 


-f(Vnf-Ua 


+ 


( 1 ) 


'Kdm = ^ J [nxdyti^ - nyd^n^ + (d^tiy - dyti^^], ( 2 ) 

where /ex is the exchange coupling and //dm is the 
Dzyaloshinskii-Moriya (DM) interaction due to the absence 
of inversion symmetry. [34-36] Here we consider that the 
DM vector D is along the bond direction in the x-y plane and 
the component along the z direction is zero, Dx = Dy = D 
and Dz = 0. [37] The field is along the z direction, = HaZ 
with a unit vector z, and the skyrmion lines are aligned in the 
same direction. In 3D systems, skyrmions are stable only in 
a small phase region in magnetic field-temperature parame¬ 
ter space [3, 38]. Here we consider a metastable skyrmion 
line lattice at zero temperature which can be obtained through 
field cooling [20] . The dynamics of spins follows the Landau- 
Lifshitz-Gilbert equation 

fiy 

dtn = ^(J • V)n - yn X Heff an x djU, (3) 
2e 

where y is the gyromagnetic ratio, a is the Gilbert damping 
constant, Heff = -^// /Sn is the effective magnetic field and J 
is the spin polarized current. The current is in the x direction 
in our simulations. 

In simulations, we calculate the skyrmion density in the 
x-y plane p(r) = n • (^j^n x dyn)/(47i), or the compo¬ 
nent. The emergent electric field induced by skyrmion mo¬ 
tion is = hn ' (Vn X dtn)/(2e). The system size is 


40 X 40 X 40 (/ex/D)^. The system is discretized into a cu¬ 
bic mesh with periodic boundary conditions in most cases. To 
count the monopole/antimonopole number, we split each sur¬ 
face of a unit cubic cell into two triangles, and calculate the 
solid angle Q/ subtended by the three spins in a triangle, 
tan(Q// 2 ) = ni • (n 2 x n 3 )[l -r ni • n 2 + 02 • 03 -r ns • ni]“^ [ 20 ] 
The monopole number of this unit cell is obtained by sum¬ 
ming the solid angle for these 12 triangles on the surface of 
the unit cubic cell, i.e. hm = (4;:)“^ Yjj=i which is quan¬ 
tized, zzm = ± 1 , 0 . 

We first study the nucleation of skyrmions in the ferromag¬ 
netic state due to the magnon instability induced by current. 
The magnon dispersion in the fully spin polarized state is 
given by [28] 

+ (4) 

2e I ^ ' 

The term J • k accounts for the magnonic Doppler shift. 
The magnon gap for / = 0 is due to the external 
field and it becomes gapless at a threshold current /^ = 
4e + 1)], signaling an instability. The averaged 

skyrmion and monopole density in the stationary state as a 
function of current is displayed in Fig. 2 (a). Right after 
the instability, skyrmion segments are nucleated accompany¬ 
ing the creation of monopoles and antimonopoles, see Fig. 
1 (c). In Fig. 2 (b), the skyrmion and monopole densities 
as a function of time are shown when the instability is trig¬ 
gered. Short skyrmion segments appear with a high monopole 
density. During the relaxation, the length of skyrmion seg¬ 
ment increases accompanying an increase in skyrmion den¬ 
sity. Once these segments meet, monopole and antimonopole 
annihilate resulting in longer skyrmion segments. This pro¬ 
cess repeats until the skyrmion lines span the whole system 
when no monopole and antimonopole is left. Meanwhile the 
skyrmion lines become straight to minimize the line tension. 
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FIG. 2. (color online) (a) Skyrmion and monopole density versus current in the stationary state for Ha - 0.8 At J - 3.84 Delfi, 

skyrmion segments are created after the instability of the ferromagnetic state, (b) Time dependence of skyrmion and monopole density when 
the instability is triggered at / = 3.1 De/h at t = 0 for Ha - 0.6 Monopoles and antimonopoles are created after the instability but they 

finally disappear in the stationary state, resulting in a straight skyrmion line lattice, (c) Dynamics of skyrmion destruction at a high current 
7 = 10.0 De/fi and Ha - 0.5 where monopoles and antimonopoles are nucleated, (d) Distribution of monopoles and antimonopoles in 

the chiral liquid phase at 7 = 8.64 De/fi and Ha - 0.8 


The monopole density decreases exponentially and the typi¬ 
cal time scale for this process is about 300 In the 

end, straight skyrmion line lattice driven by current moves as 
a whole, see Fig. I (b). 

The skyrmion density increases stepwise with current. For 
a low skyrmion density, the system is close to a ferromagnetic 
state and one would expect a similar magnonic instability as 
that in Eq. (4). For 7 close to 7^, the skyrmion density is low 
and the magnetic instability persists at the same current. In 
this case skyrmion segments do not relax into straight lines. 
When the skyrmion density becomes high, the instability cur¬ 
rent increases. The magnonic instability disappears when the 
skyrmion lattice is reached at a high skyrmion density for an 
intermediate 7. This accounts for the stepwise increase of 
skyrmion density. In these steps, no monopole exists because 
the skyrmion line lattice spans the whole system. 

At a high current, the skyrmion density decreases sharply 
indicating destruction of skyrmions [see Fig. 2 (a)]. The 
system evolves into a chiral liquid displayed in Fig. I (d), 
where spin chirality fluctuates strongly in space and time. In 
the course of skyrmion destruction, skyrmion lines break into 
segments with a monopole and an antimonopole at ends (sim¬ 
ilar to the Dirac strings). In other words, skyrmion lines are 
destroyed due to the proliferation of monopole-antimonopole 
pairs in this dynamic phase transition. The time evolution of 
monopole and skyrmion density is presented in Fig. 2 (c). 
Skyrmion density decreases while the monopole density in¬ 
creases. In the chiral liquid, both densities fluctuate justi¬ 
fying the liquid nature of the final state. These fluctuations 
are induced dynamically by current. As shown in Fig. 2 (d), 
the monopoles and antimonopoles are distributed randomly in 
space, and can be regarded as a monopole gas. The dynam¬ 
ical phase transitions associated with the skyrmion creation 
and destruction are of the flrst-order because there is a strong 
hysteresis. [28] Both transitions can be detected from current- 
voltage characteristics because the emergent electric fleld de¬ 
pends on the skyrmion density. 

We next discuss the controlled creation of monopoles and 
antimonopoles by driving part of the skyrmion lines. For 


ease of discussion, let us consider the case that current only 
flows at the bottom surface of the crystal. The skyrmion lines 
near the surface are driven by the Forentz force due to the 
current, J^(z), and start to move. This causes distortion of 
the lines and drags the skyrmion lines inside the crystal to 
move due to the elastic energy of the skyrmion lines [see 
Fig. 3 (b)]. For a small current, the elastic force can balance 
the different Forentz forces acting on the surface and mid¬ 
dle of the skyrmion line, and the skyrmion lines move as a 
whole. For an elastic skyrmion line lattice, the equation of 
motion for the lattice can be described by a displacement vec¬ 
tor u(z) = [Wjc(z), Uy{z), 0] [39, 40] 

— [tj\(z) + v(z) xz] = —£ X Jb(z) -f ^d^uiz), (5) 
y c 2 ^ 

where v = [dtUx, dtUy, 0] is the skyrmion velocity, and r] 

1 is the dimensionless viscosity due to the Gilbert damping 
for localized spins and the Ohmic dissipation by conduction 
electrons. Here i is the tilt modulus of the skyrmion lattice and 
it becomes skyrmion self-energy per unit length for a single 
line. From Eq. (5), the velocity component parallel to the 
current is V|| = -yfiJI{2e{\ -f ^^)L^] and that perpendicular to 
the current is = yfir/J/[2e(l -f ^^)L^] with the system 
size along the z direction. For a small current, the equation of 
motion is the same as that for a single skyrmion in thin Aims. 
The displacement fleld is w oc 7z^. 

However for a strong current, the skyrmion line breaks. 
The skyrmion segments at the surface and inside the crystal 
travel with different velocities. This decoupled motion of the 
skyrmion segments is possible because the monopole and an¬ 
timonopole pairs are allowed for the emergent magnetic flux 
line associated with the skyrmion. 

The above picture is confirmed by our simulations. We 
consider a T-shaped crystal shown in the inset of Fig. 3 
(a). The current is injected through the bottom of the crys¬ 
tal. The skyrmion lines produce an emergent electromagnetic 
fleld which affects the motion of conduction electrons. How¬ 
ever, as shown in the Supplement Materials [41], the emer¬ 
gent electric fleld is much smaller than the real electric fleld 
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FIG. 3. (color online) (a) Emergent electric fields at the top and bottom surfaces of the crystal. At a small current, at the opposite surfaces 
is the same indicating an elastic motion of skyrmion lines. Beyond a threshold current, at the bottom surface under a current drive is much 
bigger than the top one, implying the breaking of skyrmion lines. Inset is a sketch showing how to inject current at the bottom surface and the 
calculated current distribution [41]. (b) Snapshot of spin configuration after the breaking of skyrmion lines, (c) Monopole and antimonopole 
density versus current, (d) Distribution of monopoles and antimonopoles after the segmentation of skyrmion lines. In this calculation, we 
have used open boundary conditions in the z direction. Here Ha - 0.5 and J - 0.26 Delfi for (b) and (d). The system size is 

4 = 2L^ = 2Ly = 100 /ex/D. 


induced by the injected current, and we neglect the emergent 
electric field in the calculation of current distribution. The 
current distribution is obtained by solving = 0, subject to 
the proper boundary conditions, with V being the voltage [41]. 
As displayed in Fig. 3 (a), below a threshold current, the elec¬ 
tric field at the top and bottom are the same, indicating motion 
of skyrmion lines as a whole. Beyond the threshold current, 
the skyrmion segments at the bottom surface move at a higher 
velocity. At the same time, monopoles and antimonopoles are 
created which are confined at the bottom surface [see Figs. 
3 (c) and (d)]. At a large current, the skyrmion density near 
the bottom surface decreases due to the dynamical destruction 
of skyrmions. This causes a rapid increase of antimonopole 
density. Meanwhile the skyrmion density in the current-free 
region does not change, which results in a higher density of 
antimonopoles than that of monopoles, see Fig. 3 (c). 

The electric field along the magnetic field direction 
is particularly interesting. Let us first consider Ef due to 
the stretching of skyrmion lines by a surface current. The 
electric field is given by = v x B^/c. For a straight 
skyrmion line, E^ = 0 because = B^z only has the z 



FIG. 4. (color online) (a) Emergent electric field along the mag¬ 
netic field direction E^ when skyrmion lines are stretched by a cur¬ 
rent at the bottom surface in a T-shaped crystal in Eig. 3 (a), (b) 
Eluctuations of Ef as a function of time in the chiral liquid phase at 
J = 10.0 De/n. Here Ha = 0.5 DV/ex- 


component. When the skyrmion lines are distorted accord¬ 
ing to Eq. (5), B^ acquires an in-plane component, B^ = 
B^ld^Ux, d^Uy, 1]/ ^J(^zUx)^ + (dzUy)^ + 1 and E^ is nonzero, 
as shown in Fig. 4 (a). The induced electric field Ef OC fi 
increases nonlinearly with current in the elastic region and 
it suddenly decreases when the skyrmion line breaks, after 
which skyrmion segments become mostly aligned with the 
magnetic field. In the chiral liquid phase, Ef fluctuates around 
zero due to the fluctuations of the density of monopoles and 
antimonopoles; see Fig. 4 (b). Therefore Ef provides a clear 
transport signature of the existence of monopoles. We remark 
that the nonuniform current distribution can also produce a 
local real electric fleld E^ along the magnetic fleld direction. 
However E^ averaged over the crystal vanishes because there 
is no net current in the magnetic fleld direction. 

The existence of monopoles can also be inferred from imag¬ 
ing measurements, such as magnetic force microscope [20]. 
For instance, one may first apply a strong current to drive the 
system into the chiral liquid phase. Then the current is re¬ 
moved and the system can be trapped in a metastable state 
with monopoles. By visualizing the skyrmion density at the 
opposite surface, one can extract the information about the 
monopoles inside the crystal. 

If we start with a skyrmion line lattice in the ground state 
and drive it, the magnon instability due to the Doppler shift 
is absent. At a high current, we still have the destruction of 
skyrmion lines in nonequilibrium similar to those in Fig. 2 (a) 
via the proliferation of monopoles and antimonopoles. 

To summarize, we show that the emergent monopoles play 
an important role in the skyrmion line creation and destruc¬ 
tion process in the presence of an electric current. The current 
density required for the monopole proliferation is of the order 
of 10^^ A/m^ for typical material parameters. The monopoles 
can also be created in a controlled way by applying an inho¬ 
mogeneous current to the crystal. Especially we predict an 
induced electric fleld component along the magnetic fleld di- 
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rection, in addition to the Hall and longitudinal electric fields. 

Acknowledgments - The authors are indebted to Achim Rosch 
and Jiadong Zang for helpful discussions. The authors also 
thank Cristian D. Batista for helpful suggestions and a crit¬ 
ical reading of the manuscript. Computer resources for nu¬ 
merical calculations were supported by the Institutional Com¬ 
puting Program at LANL. This work was carried out under 
the auspices of the National Nuclear Security Administration 
of the US DOE at LANL under Contract No. DE-AC52- 
06NA25396 and was supported by the LANL LDRD-DR Pro¬ 
gram. 


* szl@lanl.gov 

[1] A. N. Bogdanov and D. A. Yablonskii, Sov. Phys. JETP 68,101 
(1989). 

[2] U. K. RoBler, A. N. Bogdanov, and C. Pfleiderer, Nature 442 , 
797 (2006). 

[3] S. Muhlbauer, B. Binz, F. Jonietz, C. Pfleiderer, A. Rosch, 
A. Neubauer, R. Georgii, and P. Boni, Science 323,915 (2009). 

[4] X. Z. Yu, Y. Onose, N. Kanazawa, J. H. Park, J. H. Han, Y. Mat- 
sui, N. Nagaosa, and Y. Tokura, Nature 465 , 901 (2010). 

[5] X. Z. Yu, N. Kanazawa, Y. Onose, K. Kimoto, W. Z. Zhang, 
S. Ishiwata, Y. Matsui, and Y. Tokura, Nature Materials 10 , 
106 (2011). 

[6] S. Seki, X. Z. Yu, S. Ishiwata, and Y. Tokura, Science 336, 198 

( 2012 ). 

[7] T. Adams, A. Chacon, M. Wagner, A. Bauer, G. Brandi, B. Ped¬ 
ersen, H. Berger, P. Lemmens, and C. Pfleiderer, Phys. Rev. 
Lett. 108 , 237204 (2012). 

[8] X. Yu, M. Mostovoy, Y. Tokunaga, W. Zhang, K. Kimoto, 
Y. Matsui, Y. Kaneko, N. Nagaosa, and Y. Tokura, Proceed¬ 
ings of the National Academy of Sciences 109 , 8856 (2012). 

[9] X. Z. Yu, Y. Tokunaga, Y. Kaneko, W. Z. Zhang, K. Kimoto, 
Y. Matsui, Y. Taguchi, and Y. Tokura, Nature Communications 
5,3198 (2014). 

[10] F. Jonietz, S. Muhlbauer, C. Pfleiderer, A. Neubauer, 
W. Miinzer, A. Bauer, T. Adams, R. Georgii, P. Boni, R. A. 
Duine, K. Everschor, M. Garst, and A. Rosch, Science 330 , 
1648 (2010). 

[11] X. Z. Yu, N. Kanazawa, W. Z. Zhang, T. Nagai, T. Kara, K. Ki¬ 
moto, Y. Matsui, Y. Onose, and Y. Tokura, Nature Communi¬ 
cations 3, 988 (2012). 

[12] T. Schulz, R. Ritz, A. Bauer, M. Haider, M. Wagner, C. Franz, 
C. Pfleiderer, K. Everschor, M. Garst, and A. Rosch, Nature 
Physics 8, 301 (2012). 

[13] L. Kong and J. Zang, Phys. Rev. Lett. Ill, 067203 (2013). 

[14] S.-Z. Lin, C. D. Batista, C. Reichhardt, and A. Saxena, Phys. 
Rev. Lett. 112 , 187203 (2014). 

[15] M. Mochizuki, X. Z. Yu, S. Seki, N. Kanazawa, W. Koshibae, 
J. Zang, M. Mostovoy, Y. Tokura, and N. Nagaosa, Nat. Mater. 
13 , 241 (2014). 

[16] J. S. White, 1. Levatic, A. A. Omrani, N. Egetenmeyer, K. Prsa, 

I. Zivkovic, J. L. Gavilano, J. Kohlbrecher, M. Bartkowiak, 
H. Berger, and H. M. Ronnow, J. Phys.: Condens. Matter 24, 
432201 (2012). 

[17] J. S. White, K. Prsa, P. Huang, A. A. Omrani, I. Zivkovic, 
M. Bartkowiak, H. Berger, A. Magrez, J. L. Gavilano, G. Nagy, 

J. Zang, and H. M. Rpnnow, Phys. Rev. Lett. 113 , 107203 
(2014). 


[18] A. Fert, V. Cros, and J. Sampaio, Nature Nanotechnology 8, 
152 (2013). 

[19] N. Nagaosa and Y. Tokura, Nature Nanotechnology 8, 899 
(2013). 

[20] P. Milde, D. Kohler, J. Seidel, L. M. Eng, A. Bauer, A. Cha¬ 
con, J. Kindervater, S. Muhlbauer, C. Pfleiderer, S. Buhrandt, 
C. Schiitte, and A. Rosch, Science 340 , 1076 (2013). 

[21] C. Schiitte and A. Rosch, Phys. Rev. B 90 , 174432 (2014). 

[22] J. Zang, M. Mostovoy, J. H. Han, and N. Nagaosa, Phys. Rev. 
Lett. 107 , 136804 (2011). 

[23] C. Castelnovo, R. Moessner, and S. L. Sondhi, Nature 451 , 42 
(2008). 

[24] M. W. Ray, E. Ruokokoski, K. Tiurev, M. Mttnen, and D. S. 
Hall, Science 348, 544 (2015). 

[25] R. Takashima and S. Fujimoto, Journal of the Physical Society 
of Japan 83, 054717(2014). 

[26] H. Watanabe and A. Vishwanath, arXiv: 1410.2213 (2014). 

[27] Y. Kawaguchi, Y. Tanaka, and N. Nagaosa, arXiv: 1508.07665 
(2015). 

[28] S.-Z. Lin, C. Reichhardt, C. D. Batista, and A. Saxena, Phys. 
Rev. Lett. 110 , 207202 (2013). 

[29] J. Fernandez-Rossier, M. Braun, A. S. Niinez, and A. H. Mac¬ 
Donald, Phys. Rev. B 69 , 174412 (2004). 

[30] V. Vlaminck and M. Bailleul, Science 322, 410 (2008). 

[31] A. Bogdanov and A. Hubert, Journal of Magnetism and Mag¬ 
netic Materials 138, 255 (1994). 

[32] J. H. Han, J. Zang, Z. Yang, J.-H. Park, and N. Nagaosa, Phys. 
Rev. B 82, 094429 (2010). 

[33] U. K. RoBler, A. A. Leonov, and A. N. Bogdanov, J. Phys.: 
Conference Series 303 , 012105 (2011). 

[34] 1. Dzyaloshinsky, J. Phys. Chem. Solids 4, 241 (1958). 

[35] T. Moriya, Phys. Rev. 120, 91 (1960). 

[36] T. Moriya, Phys. Rev. Lett. 4 , 228 (1960). 

[37] For Dx - Dy - D^- D, the ordering vector of the helical phase 
is along the [111] direction. Skyrmion lines are along the [112] 
direction when the magnetic held is along the same direction. 
The magnitude and direction of D do not change the physics 
discussed here qualitatively. 

[38] S. Buhrandt and L. Fritz, Phys. Rev. B 88, 195137 (2013). 

[39] S.-Z. Lin, C. Reichhardt, C. D. Batista, and A. Saxena, Phys. 
Rev. B 87, 214419 (2013). 

[40] J. Iwasaki, M. Mochizuki, and N. Nagaosa, Nature Communi¬ 
cations 4 , 1463 (2013). 

[41] See supplemental material at (xx) for the calculations of the 
current distribution in the crystal. 

[42] A. Neubauer, C. Pfleiderer, B. Binz, A. Rosch, R. Ritz, P. G. 
Niklowitz, and P. Boni, Phys. Rev. Lett. 102, 186602 (2009). 

[43] A. E. Petrova, V. N. Krasnorussky, A. A. Shikov, W. M. Yuhasz, 
T. A. Lograsso, J. C. Lashley, and S. M. Stishov, Phys. Rev. B 
82, 155124 (2010). 



6 



20 40 60 80 100 120 140 


FIG. 5. (color online) T-shaped crystal and the resulting current dis¬ 
tribution. The vector field denotes the x and y components of J and 
color represents the intensity. 


Calculation of the current distribution in a T-shaped crystal 


We consider a T-shaped crystal sketched in Fig. 5. A cur¬ 
rent is applied into the crystal through two electrodes at the 
bottom, and V~. The presence of skyrmion affects the 
motion of conduction electrons and changes the conductiv¬ 
ity of the crystal. The skyrmion modifies the conductivity in 
two different ways. The skyrmion produces an effective mag¬ 
netic field acting on the conduction electrons and changes the 
conductivity. Meanwhile when the skyrmions start to move, 
the motion of skyrmions produces an emergent electric field 
which also modifies the conductivity. For the transverse con¬ 
ductivity, the former contribution is known as the topological 
Hall conductivity and the latter as the skyrmion flow conduc¬ 
tivity. For MnSi, both contributions due to the skyrmions are 
found to be small compared to the background conductivity 
from the conduction electrons in experiments. [12, 42] Be¬ 
low we provide an explicit estimate of the emergent electric 
field induced by skyrmion motion and compare it to the real 
electric field. The emergent electric field is [22] 
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where as is the linear size of skyrmion and v is the skyrmion 
velocity, which can be estimated using the Thiele’s equation 
of motion for skyrmions [39] 


with a the lattice constant for the spin system. The ratio of 
to the real electric field E = pJ with p being the resistivity is 

E^ ha^ 

T - 

For MnSi, we have as ^ 10 nm, a ^ 0.45 nm and p ^ 3 x 
10“^ Q • m [43], thus we estimate E^/E 0.03. Therefore, 
we can neglect the skyrmion contribution in the calculations 
of the electric current distribution in the crystal below. 

The current distribution is governed by the equation for 
voltage V 


vV = 0. (9) 

The boundary condition in the region without electrodes is 
that the spatial derivative of the voltage normal to the bound¬ 
ary vanishes, d\V = 0. For the boundary in contact with elec¬ 
trodes, the boundary condition is V = and V = V~ =0. 
The resulting current distribution is J = -crVV, with cr be¬ 
ing the conductivity. We solve Eq. (9) numerically for the 
geometry shown in Fig. 5. 

The current mainly flows in the bottom region of the crys¬ 
tal. The current distribution does not depend on the intensity 
of the current because we can obtain the desired current inten¬ 
sity simply by rescaling . The region where current extends 
in the vertical direction is of the same order as the lateral di¬ 
mension of the crystal at the top. The nonuniform electric 
current distribution yields a local electric field in the verti¬ 
cal direction E^. However, the averaged over the crystal 
vanishes, f dPE^ = f dr^JzIcr = 0, because there is no net 
electric current in the vertical direction. 

Knowing the current distribution in the crystal, we then 
solve the Landau-Lifshitz-Gilbert equation numerically. 





































